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I. INTRODUCTION 



A. HISTORICAL 

For the past two years research has been carried out at the Naval 
Postgraduate School Department of Oceanography to study the shallow 
water interactions of various oceanographic processes that occur at 
small scales. Prior to the conduct of this experiment in May 1973, 
experiments were made in October 1971 and June 1972 from the Naval 
Undersea Center (NUC) Oceanographic Research Tower off Mission 
Beach, San Diego, California. Previous experiments included the 
simultaneous measurements of several parameters such as surface 
waves, water particle velocities, and fluctuations of temperature, 
sound speed, sound phase, and salinity. The results of these exper- 
iments are presented for the most part in MS theses of Naval Post- 
graduate School students. Smith [1971], Rautman [1971], Fitzgerald 
[1972], and Alexander [1972] reported on sound amplitude and phase 
fluctuations. Bordy [1972] made wave and water particle motion 
measurements. Krapohl [1972] made more detailed wave and turb- 
ulence measurements. Seymour [1972] initiated work on salinity 
measurements which were completed by Frigge [1973]. Duchock [1972] 
attempted to correlate sound velocity measurements with temperature 
and other fluctuations. Haley [1973] re-examined the October 1971 
experiment in an attempt to interrelate all parameters. Gossner [1973] 
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compared measured sound speed at several depths with sound speed 
computed from empirical equations, Whittemore [1973] made pro- 
gress in relating the near surface temperature field to wave motion. 

From these studies, it became apparent that a more detailed 
examination of the small scale temperature field was required before 
further progress could be made in determining the interactions of 
other parameters. 

B. STATEMENT OF THE PROBLEM 

The creation of the small scale temperature field in shallow water 
is a complex process and may occur at various depths in the water 
column. At the surface, wind stresses create waves and if the stress 
is sufficient to produce wave breaking (white caps), then a turbulent 
transfer of momentum downward will result with an attendant down- 
ward flux of any excess heat in the surface water. Wind stress will 
also produce surface drifts. The shears from these drifts produce a 
certain amount of turbulence which may also be transmitted downward. 
This process of surface originated turbulent mixing extends downward 
to a depth proportional to the wavelength of the surface waves. 

The ocean bottom also exerts a stress on water moving over it 
due to, for instance, tidal currents. The resulting shear zone causes 
turbulent mixing which affects the small scale temperature field and 
extends upward a distance proportional to the stress exerted. 
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The presence of a thermocline inhibits the vertical extent of the 
turbulence caused by the processes described above. Since the degree 
of stability of the water column tends to define the degree of vertical 
motion, the increase in stability implied by the density gradient will 
reduce the vertical motion of the turbulence and thus the vertical 
mixing. 

It is paradoxical that the thermocline, which inhibits the turbulent 
mixing due to surface and bottom shears, may itself be the point of 
creation of turbulence. Internal waves propagate along the thermocline 
and can lead to the creation of intense turbulence. If the water masses 
separated by the thermocline experience relative motion, then the shear 
across the thermocline will produce turbulence. This relative motion 
may be caused by tidal currents or wind drifts. [LaFond, 1962] 

It is interesting that turbulent mixing from surface and bottom 
shears tends to intensify the thermocline whereas turbulence caused 
by internal waves or current shears tends to dissipate it. Wind 
induced mixing redistributes heat downward to the thermocline so 
the deep water becomes warmer; bottom mixing essentially re- 
distributes heat downward away from the thermocline so the upper 
level becomes cooler. The gradient thus becomes stronger. The 
internal waves tend to weaken the gradient through complex non-linear 
wave interactions and make the temperature difference smaller. 

The greater the shear across the thermocline, the greater the 
possibility of instability and turbulence. It is well known that the 
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intensity of the shear determines the rate of growth of an unstable 
disturbance in a steady, stable, parallel, shearing motion. The shear 
instability may be described by the local Richardson number, J(z), 




where g is the acceleration due to gravity, £ is the reference 
density, Mt is the rate of density change across the thermo- 

cline, and is the velocity shear across the thermocline. When 

the Richardson number is less than a critical value (J(z)^ 0. 25) then 
turbulence is possible though not certain [Phillips, 1969]. Since the 
maximum shear is normally found at the crests and troughs of waves, 
this turbulence should be expected to occur most of the time or with 
greatest intensity at these points. If the thermocline is intense, the 
critical shear will be reached earlier than if it were diffuse; an 
internal wave in an intense thermocline will dissipate its energy into 
turbulence more quickly than one in a diffuse thermocline. The tur- 
bulent mixing resulting from the passage of an internal wave results 
in the reduction of the shear in the thermocline to a value less than 
the critical value. Creation of turbulence by this means will then 
presumably cease until the thermocline is again restored. 

An internal wave need not be a turbulent phenomenon. If there is 
no turbulence being created at the thermocline, then the regenerating 
forces of turbulent mixing from surface and bottom shears may 
restore the intensity again. 
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The interplay between propagation of energy in an internal wave 
and its dissipation may be summarized in the following way. An 
internal wave is initiated by some process and propagates along a 
density gradient. The shear at the boundary increases possibly due 
to the wave entering shallow water and changing its form or possibly 
due to an increase in the intensity of the density gradient or thermo- 
cline. If this shear exceeds a critical value, then turbulence may be 
created, probably at the crest or trough. The resulting turbulent 
mixing reduces the gradient possibly below the critical value. The 
wave then either continues propagating with energy reduced by losses 
to turbulence, or, if the shear is not reduced below the critical value, 
the wave is entirely dissipated. Areas of density gradients reduced 
by the mixing of past waves will not be able to support waves of such 
intensity of shear. 

Theoretically, then, there is a cycle between an intense thermo- 
cline with intense internal wave activity producing turbulence, and a 
diffuse thermocline with internal waves producing little or no 
turbulence. 

The kinetic energy of the turbulence is eventually dissipated by 
entrainment. Shears at the boundaries of the turbulent eddies formed 
by the mixing processes described above cause the shedding of smaller 
eddies, among which the energy of the original eddy is distributed. 

The energy is cascaded through eddies of decreasing size and energy 
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until eventually the energy is no longer sufficient to overcome the 
viscous forces of the surrounding water. The remaining energy is 
then dissipated as heat energy. 

Small scale temperature gradients may exist in a field even 
though all velocity gradients have been removed by viscosity because 
the scales at which molecular diffusion of heat becomes effective are 
much smaller than those at which molecular viscosity is effective. 
This small scale temperature field is subject to being moved about 
as a passive scalar riding on the water particle motion caused by 
surface and internal waves. 

C. OBJECTIVES 

An experiment was conducted in May 1973 in order to make a 
careful examination of the small scale temperature field in the near 
surface region where the influence of surface waves is felt. The 
principal sensors were a moveable array of thermistors, a wave 
gauge, and a flowmeter. The object was to verify Whittemore's 
suggested wave- temperature coupling mechanism. 

Particular attention was given to determining time and length 
scales. The scales that can be resolved are determined by the char- 
acteristics of the sensors and their spatial relationships with each 
other. For this experiment spatial scales were determined by the 
minimum thermistor spacing (about eight centimeters) and by the 
maximum spacing (about two meters). Minimum temporal scales 
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depended on the response time of the thermistors (0. 015 seconds). 
Maximum scales depended on the length of the record. A maximum 
scale of five percent of the record length of about 20 minutes was used 
for the maximum lag time in calculating correlation functions giving 
an upper limit of one minute period waves. 

The conditions of the near surface region were fulfilled by the 
shallow water depth of approximately 19 meters. 

The May date was chosen partly for logistic reasons but also 
because only moderate wave activity was expected and vertical tem- 
perature gradients would not be intense. Internal wave activity was 
not expected to be a major source of turbulence. Unfortunately, 
conditions did not occur as predicted. Wind waves and swell were 
small or nonexistent and internal waves dominated the temperature 
scale much of the time. Although these conditions detracted from the 
original aim of the experiment, they provided an excellent opportunity 
to examine small scale (turbulent) temperature fluctuations in the 
presence of internal waves. 

This thesis presents an analysis and discussion of the data taken 
during the May 1973 experiment conducted from the NUC tower. 
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II. DESCRIPTION OF THE EXPERIMENT 



A. LOCATION 

The experiment was conducted from the Naval Oceanographic 
Research Tower operated by the Naval Undersea Center (NUC) 
located about one mile off Mission Beach, San Diego, California. 

The tower consists of four deck levels about 40 ft square each on a 
frame of four legs. One of the decks contains a small deck house 
suitable for housing electronic equipment. Vertically mounted rails 
from one of the lower decks run down three sides of the support frame; 
equipment may be mounted on these rails and lowered to any pre- 
selected depth. The rails on the seaward face were used for the 
experiment since wave action was expected to come from that 
direction. The tower offers stable support and ample working space. 
Water depth at the tower is about 18. 3 m with a smooth sandy bottom. 
[La Fond, 1965] 

B. TIME 

The experiment was conducted 15 and 16 May 1973; this time of 
year is a time of transition between winter and summer conditions. 

Air and sea surface temperatures are normally in the vicinity of 
16.7°C, low wind, sea, and swell activity is usual, the thermocline 
is not likely to be strong, and internal wave activity is not yet intense. 
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Actual meteorological and oceanographic environmental factors 



observed were: 



Day 1 : 



RH: 

SST: 



air temp: 
dew point: 



o 

15.5 C 
13. 9°C 
84% 0 



16.5 C 



overcast 
wind: calm 

sea sfc: negligible wind 

waves; swell 0. 7 m from west 



Day 2: 



RH: 

SST: 



air temp: 
dew point: 



14. 4- 15 . 5°C 
13. 3-13. 9°C 
-84% 




overcast 
wind: calm 

sea sfc: capillary waves; 

swell 0.6 m from west 



Conditions were essentially constant over the two days except the sky 
cleared at the end of the last run on the second day. Nansen casts 
were made for all runs, results are in table I; bathythermograph and 
sea surface temperature measurements were made for each run with a 
shallow BT and a bucket thermometer. The calibration of the BT was 
uncertain and there was considerable hysteresis. However, an 
estimation of thermal gradients was made for each run; the results 
are in table II and figure 1. 

Tidal conditions were as follows: 



Day 1 



1340 low 
1955 high 



1125 ebb 
1729 flood 



1407 slack 



Day 2 



0917 high 
1409 low 



0933 slack 
1432 slack 



1158 ebb 
1754 flood 
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TABLE I 



RESULTS OF NANSEN CASTS 



RUN 








NUMBER 


TIME 


DEPTH 


SALINITY 


1 


1315 


3 m 


33. 743 0/00 


2 


1440 


8 m 


33. 736 0/00 


3 


1600 


14 m 


33. 739 0/00 


4 


1730 


3 m 


33. 738 0/00 


6 


1100 


3 m 


33. 743 0/00 






TABLE II 






ESTIMATION 


OF THERMAL 


GRADIENTS 


RUN 




TOP OF 


INTENSITY OF 


NUMBER 




GRADIENT 


GRADIENT 


1 




4. 25 m 


0.98 C/m 


2 




10. 5 m 


0. 88° C/m 


3 




7. 6 m 


0.79° C/m 


4 




14. 7 m 


0.25° C/m 


6 




sfc 


0.44° C/m 
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Temperature (°C) 




Temperature (°C) 




( I indicates depth of frame) 



Figure 1. Bathythermograph traces 
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Depth(m) Depth (m ) 

18.3 0 18.3 



C. EQUIPMENT USED 



A specially constructed 6x6x8 foot frame was mounted on the 

vertical rails of the seaward face of the tower by means of sliding 

shoes. An array of seven individually movable thermistors and a two 

channel electromagnetic flowmeter were mounted on the outside face 

of the frame. The frame was then lowered to the predetermined depth 

by means of an electric winch mounted on the deck at the level of the 

railhead. Approximate depth was determined by means of a marked 

cord attached to the top of the frame and payed out from the winch 

deck. An I-beam was mounted horizontally on the seaward face of 

the deck house extending over the equipment frame. A trolley ran on 

the underside of the beam; a Baylor wave gauge was suspended between 

the trolley and the top of the outside face of the equipment frame 

[See Fig. 2]. Since the frame was displaced horizontally as it slid 

o 

down the rails due to a 5 cant in the rails, the trolley was moved to 
keep the wave gauge vertical and directly over the flowmeter and ther- 
mistor array. Cable runs ran from all equipments to the recording 
room in the deck house. 

The core of the experiment was seven thermistors. The ther- 
mistors were type K496 from Fenwall Electronics Inc. , designed for 
oceanographic work. The thermistors consist of two specially aged 
beads in series; since there is a statistically small probability of both 
beads drifting or drifting in the same direction with time, they are 
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exceptionally stable and have a high repeatable accuracy. The thin 
glass watertight sealant around the beads allows a short response 
factor of 0.015 sec. High resolution is obtained from the high resist- 
ance (15 K ^ at 25°C) and a high resistance to temperature change 
o o 

ratio (540 K ft / C at 25 C). Due to low self heating, measurements 
accurate to the nearest 0. 01°C can be obtained in seawater [Shonting, 
1968]. The thermistors were mounted in a 0. 3x9 cm aluminum tube 
with the beads extending just beyond the lip of the tube. The tube was 
then set in a 6 cm diameter epoxy base in which was sealed the con- 
nection between the cable run and the thermistor leads. The beads 
were protected in transit by a stainless steel sleeve that fitted over the 
tube; during the experiment they were protected from fish and seaweed 
by a 3 mm diameter bronze cage (Plate 1). Pipe straps were attached 
to the epoxy bases which allowed the thermistors to be mounted on the 
frame. 

Associated with each thermistor was a Wheatstone bridge. Of the 
four legs of the bridge, two were occupied by equal fixed resistors; 
the balancing leg consisted of a 1 turn and a 10 turn potentiometer 
for coarse and fine tuning; the thermistor formed the fourth leg 
(Fig. 3). The bridge was constructed so that it could be balanced 
with the thermistor in place. A decade box could then be substituted 
for the thermistor and the bridge rebalanced by means of the decade 
box. The balancing resistance and thence the temperature of the 
thermistor would then be known. With the thermistor back in place, 
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Plate 1 

Thermistor and Guard 
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Battery, 6.5 V 




Figure 3. Circuit Diagram for Wheatstone Bridges 
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Calibrate 



any voltage fluctuation across the bridge would be a function of the 
thermistor temperature. Power was supplied by a 6. 5 V mercury 
battery. 

The thermistors were calibrated by immersion in a fresh water 

bath of known temperature. The temperature of the bath was measured 

o 

by a quartz thermometer which had itself been calibrated at 0 C in a 

distilled water and ice bath. Plots were made of bridge output versus 

o o 

temperature for temperatures between 10 C and 21 C in one degree 
increments. The curves were sufficiently linear for a valid assump- 
tion of linearity over a range of a few degrees (Fig. 4). 

The calibration factor for the individual thermistors was then 
determined by solution of the Wheatstone bridge equations: 



E = (R D +R 1 )i 2 



E— ( Rt^"R2 H j 



V 1 2 R £) - 1 i R t 



or 

y_ r D E _ r T E 

(Rq+Rj R,p+R 2 

where R^ is the instantaneous resistance of the thermistor, R^ is the 
balance point of the bridge as determined by the decade box, R^ and 
R^ are the constant resistances, and E is the applied voltage, and V 
is the output of the bridge when unbalanced. The calibration factor 
thus determined was typically 14. 9°C/V, (Fig. 5). 
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Resistance (KJl) 




Temperature (°C) 

Figure 4. Typical Temperature vs 
Resistance Calibration Curve for Thermistors 
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Figure 5* Calibration Factor for Wheatstone Bridge 
Balanced at 14°C 
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The Baylor wave gauge used has an accuracy of within 0.01 of 
wave height and a response time of less than 0. 06 sec. The gauge 
used was 4. 1 m long; its calibration factor was 7. 7 m/V. The cal- 
ibration curve is figure 6. 

The two channel electromagnetic flowmeter was an Engineering 
Physics Company Water Current Meter model EMCM3B. The ac- 
curacy is to within 0. 015 m/sec and has a flat response out to at 
least 0.5 seconds. The calibration factor was 0.541 m/sec/V. A 
complete description of this flow meter and calibration is given by 
Steer [1972], 

D. RECORDING 

All signals except those from the thermistors were recorded 
directly from the instruments. The Wheatstone bridges had to be 
balanced with the thermistors in the ocean during the procedures 
described for calibration. Voltage fluctuations caused by the tem- 
perature fluctuations in the ocean were never serious enough to 
prevent balancing the bridges. 

Since the thermistor signals were very small (on the order of 
millivolts), they had to be amplified before being recorded. Seven 
identical amplifiers were not available, so three Preston 8300XWB 
and four Hewlett-Packard 2470B amplifiers were used. Unfortunately, 
the only gain setting they had in common was xlO. For some runs in 
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Figure 6* 



Calibration Curve for Baylor Wave Gauge 



water with little thermal activity a higher gain setting could be used: 
x30 on the Preston amplifiers and x50 on the Hewlett-Packards. In 
the presence of internal waves when temperature excursions were 
relatively great, the lower gain setting was used to avoid overdriving 
the tape recorder. 

The Preston amplifiers were equipped with a filter; consequently, 
on the first day channels with the Preston amplifiers were filtered at 
10 Hz low pass to filter out any 60 Hz noise. However, on the second 
day, the filters were not used in order not to phase shift the signals. 

All signals were then recorded simultaneously on a Sangamo 
3562 14 channel FM magnetic tape recorder. Full scale deflection 
was set at + 1 V rms. Tape speed used was 17/8 ips; this allowed 
maximum speedup (to 60 ips) during playback. 

E. EXPERIMENTAL PROCEDURE 

Equipment set up for the two days was similar in most respects. 

On the first day the thermistor array was vertical as shown in figure 7. 
On the second day the array was horizontal for the first two runs 
(Fig. 8), then for the last three runs, thermistor one was moved 
below the remainder of the horizontal array to give some indication 
of vertical temperature gradient. 
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o 
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Thermistors 






*7 
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Figure 7. Configuration of Array for 
Runs One, Two, Three, and Four 
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Figure 8. Configuration of Array for Run Six 
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162.5 



Run times and depth of the top of the frame are given below: 



Run number 

1 

2 

3 

4 

5 

6 

7 

8 
9 



Time 

1300-1340 

1437-1510 

1607-1640 

1722-1755 

1017-1034 

1050-1201 

1244-1315 

1339-1400 

1410-1430 



Depth (m) 

2.7 
5.2 
11 . 3 

2.7 

2.7 

2.7 
8.5 
13 . 1 
13 . 1 
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III. DATA ANALYSIS 



A. DATA PREPARATION 

Analysis of the data was done at the Naval Postgraduate School. 
Initially all signals were graphically displayed on a strip chart for 
visual inspection, polarity checks, and checks for errors or anomalies. 
Immediately, two runs were discarded. Runs five and eight had been 
started with amplifier gains too high and too low and had overdriven 

I 

and underdriven the tape recorder respectively. Runs six and nine 
had been made at the same depths and with the same configuration and 
therefore replaced five and eight. On runs seven and nine, there was 
a great deal of noise on the signal from the sixth thermistor; these 
runs were not analyzed. 

Some long term fluctuations were noticed on all runs. Two runs 
showed significant internal wave activity. Run three showed a distinct 
internal wave in the top three thermistor traces (Pig. 9). Run six 
was started in the middle of an internal wave, spanned a second distinct 
wave, and included two patches of increased turbulence that may have 
been the remains of other internal waves (Figs, 10, 11). 

Twenty minute portions of each run were chosen for analysis. 

This was the first 20 minutes in all but run six. In run six the portion 
for analysis did not include the large internal waves; the microstructure 
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in the internal waves did not appear to be representative of the rest 
of the run. 

Once suitable portions of the data had been identified, they were 
digitized on a hybrid computer CI5000, XDS9300. Some runs required 
filtering of 60 Hz noise; this was done with Krohn-hite 3750 filters. 
Playback was at 60 ips with a sampling rate of 400 Hz, corresponding 
to a real time sampling rate of 12.5 Hz. The signals were further 
amplified during digitization to take advantage of the full dynamic 
range of the computer and thus increase resolution. Amplification was 
either x50 or xlOO depending on the amount of fluctuation in the signal. 

The seven track digital tape thus made on the hybrid computer 
was then transformed onto a nine track digital tape compatible with 
the school's IBM 360-67 computer on which the analysis was 
accomplished. 

Since the scales to be analyzed were on the order of about one to 
100 seconds, it was decided to filter out digitally any activity with 
a period of 100 seconds or greater. A Fourier analysis was made 
of each 20 minute segment of data; the coefficients were used to 
generate a time series describing the waves of period 100 seconds and 
greater. This series was then subtracted from the original time series. 
It was the filtered time series that was used for analysis. 

The DC component and any trend were then deleted. 
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B. METHOD OF ANALYSIS 



Temporal auto- and cross-correlations were computed for selected 
sensors. A Parzen window was first applied to the correlation func- 
tions. Auto- and cross-spectra were computed by a Fourier trans- 
formation on the correlation functions. Coherence and phase were 
computed from the cross-spectra. Due to the desire to cut down on 
computation time but keep the 20 minute record length, the number of 
data points was arbitrarily cut in half by deleting every other data 
point. This reduced the sampling rate to 6.25 Hz resulting in a 
Nyquist frequency of 3. 1 Hz. The reduction in the Nyquist frequency 
was not considered critical, since the spectra became essentially 
flat and the signal lost in noise before the 3. 1 Hz frequency was 
reached. Aliasing was experienced possibly from 60 Hz noise which 
was never effectively filtered out. The energy was always easy to 
identify, however, and appeared in narrow spikes at 2. 59 Hz, and 
occasionally at 2. 09, 1. 04, and 0.67 Hz. 

Spatial correlation was computed by averaging cross variance 

(C. V. ) between the first thermistor (Th ) and the others (Th ), and 

I i 

normalizing with the variance of the first thermistor, 

c v = 1/T / Thl(t,x)Th-(t, x+ A Xj) dt i=1 _ 7 

1/T J (Th^t/x)) 2 dt 

where T is the time over which the cross variance was computed 
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(20 min). The spatial correlation was then plotted against separation 
distance between thermistors using thermistor one as a reference; 
curve fitting was done by linear interpolation. The unnormalized 
spatial correlation functions were then transformed to obtain wave 
number power spectra in the same manner as described for the 
temporal functions. 



i 



40 



IV. RESULTS OF THE ANALYSIS 



A. TEMPORAL SCALES 

One of the most common ways of presenting scales of inhomo- 
geneities is to scale them by correlation times. In the present case, 
the correlation time is defined as the time required for the normalized 
autocorrelation function to fall to a value of 1/e (=0. 368). In water 
with a constant temperature gradient (such as the mixed layer), it 
might be expected that the shortest correlation times would be at the 
surface where turbulence and short period waves dominate the structure. 
With increasing depth where turbulent mixing is not so intense and 
small scale temperature gradients have begun to smooth out, the 
correlation times should lengthen out to reflect the dominance of the 
surface wave motion. In the vicinity of a thermocline having internal 
waves, the times should again increase to reflect the dominance of the 
longer period internal waves, unless the internal waves are creating 
turbulence, in which case the correlation times can shorten again. 

In the shallow water of this experiment, turbulent mixing (and there- 
fore shortening of the correlation times) may also have been caused 
by tidal currents. Consequently it is not surprising that the temporal 
correlations shown in table III show no dependence with depth. 

In addition, the vertical dimensions of the thermistor array were 
small with respect to the extent of the water masses involved. The 
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